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Abstract. The recent detection of gamma-ray lines from radioactive ^^Al and ®°Fe in the Milky Way by the 
RHESSI satellite calls for a reassessment of the production sites of those nuclides. The observed gamma-ray line 
flux ratio is in agreement with calculations of nucleosynthesis in massive stars, exploding as SNII (Woosley and 
Weaver 1995); in the light of those results, this observation would suggest then that SNII are the major sources 
of 2'^Al in the Milky Way, since no other conceivable source produces substantial amounts of ^''Fe. However, more 
recent theoretical studies find that SNII produce much higher ^'°Fe/'^^Al ratios than previously thought and, 
therefore, they cannot be the major ^® Al sources in the Galaxy (otherwise ""Fe would be detected long ago, with 
a line flux similar to the one of ^^Al). Wolf-Rayet stars, ejecting ^^Al (but not ^°Fe) in their stellar winds, appear 
then as a most natural candidate. We point out, however, that this scenario faces also an important difficulty. 
Forthcoming results of ESA's INTEGRAL satellite, as well as consistent calculations of nucleosynthesis in massive 
stars (including stars of initial masses as high as 100 Mq and metallicities up to 3 Zq), are required to settle the 
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1. Introduction 

^^Al is the first radioactive nucleus ever detected in the 
Galaxy through its characteristic gamma-ray line signa- 
ture, at 1.8 MeV (Mahoney et al. 1982). Taking into 
account its short lifetime (~1 Myr), its detection con- 
vincingly demonstrates that nucleosynthesis is still ac- 
tive in the Milky Way (Clayton 1984). The detected flux 
(~4 10~^ cm~^ s~^) corresponds to ~2 M0 of ^^Al cur- 
rently present in the ISM (and produced per Myr, as- 
suming a steady state situation). The COMPTEL instru- 
ment aboard CGRO mapped the 1.8 McV emission in 
the Milky Way and found it to be irregular, with promi- 
nent " hot-spots" probably associated with the spiral arms 
(Diehl et al. 1995). The spatial distribution of ^""Al sug- 
gests that massive stars are at its origin (Prantzos 1991, 
1993, Prantzos and Diehl 1996). However, it is not yet 
clear whether the majority of observed ^^Al originates 
from the winds of the most massive stars (i.e. above 30 
Mq, evolving as Wolf-Rayet stars) or from the explosions 
of less massive stars (i.e. in the 12-30 Mq range, explod- 
ing as SNII) ; the uncertainties in the corresponding stellar 
yields are still quite large (see Sec. 2) and do not allow to 
conclude yet. 

Clayton (1982) pointed out that SNII explosions pro- 
duce another relatively short lived radioactivity, ®"Fe (life- 



time ^2 Myr). Since WR winds do not eject that iso- 
tope, the detection of its characteristic gamma-ray lines 
^ in the Milky Way would constitute a strong argument 
for SNII being at the origin of ^^Al. Based on detailed 
nucleosynthesis calculations of SNII (from Woosley and 
Weaver 1995) Timmes et al. (1995) found that the ex- 
pected gamma-ray line flux ratio of ^°Fe/^^Al (for each of 
the two lines of ^"Fc) is 0.16, if SNII arc the only sources 
of 26 Al in the Milky Way. 

The Reuvcn Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) detected the galactic ^^Al emission at 
a flux level compatible with previous observations (Smith 
2003a). Most recently. Smith (2003b) reported the first 
ever detection of the Galactic ^''Fe gamma-ray lines with 
RHESSI; their combined fluxes correspond to a signifi- 
cance level slightly higher than 3 tr. The line fiux ratio 
6°Fe/26Al is found to be 0.16 (for each ^°Fe fine), precisely 
at the level predicted by Times et al (1995) on the basis of 
Woosley and Weaver (1995) nucleosynthesis calculations. 

This finding of RHESSI appears as an impressive con- 
firmation of a theoretical prediction. However, more recent 
studies of SNII nucleosynthesis produce different values for 
the ^^Fe/^^Al ratio (see next section), considerably higher 
than the one of Timmes et al. (1995). Combined with the 
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^ At 1.117 and 1.332 MeV, resulting from the decay of its 
daughter nucleus ^°Co 
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RHESSI finding, the new theoretical results call for a re- 
assessment of the ^^Al sources in the Milky Way. In this 
work we discuss those results and their implications. We 
argue that none of the proposed sources of ^^Al satisfies 
all observational constraints at present. Forthcoming ob- 
servations by the INTEGRAL satellite, combined with a 
new generation of stellar nucleosynthesis models (for ro- 
tating massive stars up to 100 Mq and metalhcities up to 
3 Z0) will probably be required to settle the issue. 

2. ^^Al and ^°Fe: revised yields and 7-ray line 
fluxes 

Four different groups (to our knowledge) have performed 
calculations of nucleosynthesis in massive stars, estimat- 
ing the amounts of both ^^Al and ^"Fe and covering a 
relatively extended grid of stellar masses: Thielemann et 
al. (1995), Woosley and Weaver (1995), Rauscher et al. 
(2002) and Limongi and Chieffi (2003). In the first case, 
however, presupernova calculations are made in pure He- 
cores and the amounts of hydrostatically produced ^^Al 
are seriously underestimated; therefore, those results are 
not discissed in the following. 

The calculations of Woosley and Weaver (1995, here- 
after WW95) and of Rauscher et al. (2002, hereafter 
RHHW02) are made with essentially the same stellar evo- 
lution code, but the latter benefit from improved stellar 
physics and, especially, an updated library of nuclear re- 
action rates. Thus, the RHHW02 results supersede those 
of WW95, at least for solar mctalhcity stars (WW95 is 
the only published work providing yields of radioactive 
nuclei for an extended grid of stellar metallicities) . Both 
those calculations take into account neutrino-induced nu- 
cleosynthesis during the supernova explosion, which in- 
creases the ^^Al yield by about 40% on average (WW95). 

Finally, the Limongi and Chieffi (2003, hereafter LC03) 
calculations are done with a different stellar evolution 
code but with essentially the same set of nuclear reaction 
rates as RHHW02 (the REACLIB library of Rauscher and 
Thielemann). They adopt a different treatment for the 
study of the explosion than RHHW02 and they do not 
take into account neutrino-induced nucleosynthesis. 

The situation concerning the ^^Al and ^"Fe yields of 
those calculations is summarized in the first and second 
panel of Fig. 1, respectively. In the top panel, it is clearly 
seen that the ^^Al yields of RHHW02 are substantially 
smaller than those of WW95, by a factor two on aver- 
age. That difference is obviously due to the different input 
physics adopted in the two studies. 

The LC03 yields of ■^^Al are even smaller than those 
of RHHW02, and that difference can be attributed, at 
least partially, to the neglect of the neutrino-induced nu- 
cleosynthesis in the former study. Note that, in order to 
account for the uncertainties of the supernova explosion 
LC03 study a range of explosion energies, and this affects 
(slightly) the ^^Al yield of their lowest mass stars. Note 
also the interesting "convergence" of the three calculations 
in the case of the 20 Mq star, perhaps because the proper- 
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Fig. 1. Yields (in Mq) of ^^Al (top) and ^°Fc{middle) 
and the ^°Fe/^^Al gamma-ray line flux ratio {bottom) as 
a function of stellar mass, according to calculations by 
WW95 {open squares connected by solid lines), RHHW02 
{filled circles connected by dotted lines) and LC03 {filled 
pentagons); in the latter case, different points for a given 
stellar mass correspond to different explosion energies. All 
calculations are for solar metallicity stars. The dotted line 
in the lowest panel corresponds to the RHESSI value of 
0.16 (Smith 2003b). 



ties of that particular stellar mass are better constrained 
after the extensive study of SN1987A. 

In the case of ^°Fe, RHHW02 obtain yields twice as 
large as WW95, on average. The reason for that discrep- 
ancy is probably the improved library of nuclear reaction 
rates of RHHW02. Combined with the results for ^^Al, 
it becomes obvious that RHHW02 get ^"Fe/^^A\ ratios 
four times larger than WW95. The corresponding results 
of LC03 are in excellent agreement with WW95 above 20 
M0 and in fair agreement with those of RHHW02. In the 
13-15 M0 range, the ^"Fc yields of LC03 depend strongly 
on the explosion energy, with the lower energies leading to 
higher yields. Note, however, that in all cases ^"Fe is pro- 
duced by successive neutron captures on Fe-peak nuclei; 
the last step (^^Fe(n,7)^''Fe involves the unstable nucleus 
(^^Fe, for which there are no experimental data concerning 
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its neutron capture cross-section. The nuclear uncertain- 
ties on its yield are thus quite important.^ 

The corresponding ratio of ^°Fe/^^Al by number (i.e. 
the yield ratio divided by 60/26) for each stellar mass ap- 
pears in the bottom panel of Fig. 1. The results of WW95 
are, on average, close to the value of 0.16 (dotted hori- 
zontal line), mentioned in the RHESSI discovery report 
of ™Fe (Smith 2003b), while those of RHHW02 and LC03 
are substantially above that value for almost all the stellar 
masses. 

To compare properly with observations, these yields 
should be convolved with a stellar Initial Mass Function 
(IMF) and we adopt here the Salpeter IMF, a power-law'^ 
with slope x=-1.35, in order to obtain the number ratio 
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as a function of the upper limit of integration Mjjp. The 
results are plotted on Fig. 2. In the case of LC03 two 
curves are shown: LC03II corresponds to the high ^°Fe 
yields (low explosion energies) and LC03L to the low ^°Fe 
yields (high explosion energies). In all cases the thick por- 
tions of the curves correspond to the stellar mass range 
covered by each study, while the filled hexagons mark the 
highest mass of each calculation and thus provide the IMF 
weighted value over the whole mass range covered by each 
study. 

The WW95 yields, integrated up to the highest mass 
of that study (M(7P=40 M©) lead to a number ratio 
™Fe/^^Al=0.18, i.e. very close to the value 0.16 advanced 
by Timmes et al. (1995) on the basis of those same yields 
and the same IMF. It is precisely that theoretical predic- 
tion, well within reach of modern instruments, that made 
^°Fe a prime target for astrophysical gamma-ray spec- 
troscopy. The RHESSI discovery apparently confirms that 
prediction. However, an inspection of the more recent re- 
sults shows that the modern theoretical expectations are, 
in fact, much higher: near unity for RHHW02 and LC03H 
and above 0.4 in the case of LC03L. In the light of those 
results, the RHESSI discovery at the level predicted by 
Timmes et al. (1995) looks more as a coincidence. 

Assuming that the recent theoretical results arc not to 
be substantially revised in the future and that the RHESSI 
result is confirmed, what are the implications for our un- 
derstanding of the origin of ^^Al? The obvious conclusion 
is that the bulk of galactic ^^Al, detected by various instru- 
ments including RHESSI, is not produced by the source 
of ^°Fe: if this were the case, then ^"Fe would be detected 
with a line flux similar to the one of ^^Al. Obviously, an- 
other source of ^^Al is required, producing much smaller 
60pg/26Al ratios than the SNII. 



The WW95 yields are calculated with cross section value 
of 1.8 mb for the neutron capture on ^^Fe, where the RHHW02 
and LC03 calculations adopt a value of 3.4 mb, see Woosley et 
al. (2003). 

^ A power-law IMF is defined (as a function of stellar mass 
M) as : $(M) = dN/dM = AM-'^^+''\ 
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Fig. 2. The expected ratio of ^°Fc/^^Al decays (for each 
of the two ^"Fe lines), convolved with a Salpeter stellar 
Initial Mass Function, is shown as a function of the up- 
per stellar mass limit of the convolution integral (see Equ. 
1). The four curves correspond to the four different sets 
of stellar yields, with their thick portions corresponding 
to the mass range covered in those works (sec text). The 
dotted horizontal line at 0.16 is the ^"Fe/^^Al ratio re- 
ported by RHESSI. Filled pentagons mark the upper mass 
in each of the four studies; all recent calculations predict 
much higher values than the older calculations of WW95 
or the observed value of RHESSI. Only by taking into ac- 
count the ■^^Al yields of massive Wolf-Rayct stars [thin 
portion of the curves beyond the masses indicated by the 
filled pentagons, obtained by adding data for WR stars 
from Meynet et al. 1997) one may obtain ^^Ye/'^^PA ratios 
compatible with the RHESSI results. 



The obvious candidate source is Wolf-Rayet stars, as 
has been argued in many places over the years (e.g. 
Dearborn and Blake 1985, Prantzos and Casse 1986, 
Prantzos 1991 and 1993, Prantzos and Dichl 1996, Meynet 
et al. 1997, Knodlscder 1999). The winds of those massive, 
mass losing stars, eject large amounts of ^^Al produced 
through H-burning in the former convcctive core, before its 
radioactive decay (in stars with no mass loss, those quan- 
tities of ^®A1 decay inside the stellar core before the final 
explosion and never get out of the star). Note that WR 
stars eject negligible amounts of ^°Fe, since that nucleus is 
produced at more advanced stages of the stellar evolution 
than ^^Al and there is no time for it to be ejected before 
the final explosion (e.g. Prantzos et al. 1987). However, no 
complete calculations of WR stars (i.e. of massive stars, 
say above 40 Mq , with mass loss and up to the final explo- 
sion) are available up to now. The calculations of Woosley 
et al. (1995) concern only the advanced evolution of mas- 
sive He cores and ignore any contribution of the WR winds 
to the ^^Al yields (besides, it is difficult to link the mass 
of their calculated He cores to the mass of the correspond- 
ing main sequence stars). Thus, the total ^^Al and ^°Fe 
yields of those stars (i.e. the sum of the masses ejected by 
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Fig. 3. Yields of SNII from WW95 for ^^Al {top) and 
^°Fc (bottom), for three different values of the initial stellar 
metallicity (as indicated in the bottom panel). 



the winds and by the explosion) is unknown at present. 
Although the number of such stars in a normal IMF is 
small, the amounts of ^^Al in the winds are extremely 
large and affect considerably the overall budget. In what 
follows, we assume that those stars eject negligible total 
amounts of ^°Fe (otherwise one faces the same problem 
with a high ^"Fe/^^Al ratio as before). 

Adopting the ^^Al yields of non-rotating WR stars 
of solar initial metallicity by Meynet et al (1997), which 
concern stars of solar metallicity in the 25-120 Mq mass 
range, and combining them with the aforementioned SNII 
yields, one obtains the ^°Fe/^^Al ratio expected by the 
total mass range of massive stars, during all the stages 
of their evolution; this is expressed in Fig. 2 by the con- 
tinuation of the four theoretical curves above the masses 
indicated by the filled pentagons. It can be seen that the 
RHESSI result is recovered in that case, provided that at 
least half of ^^Al originates from WR stars (in the case 
of LC03L), or even that 80% of ^^Al originates from WR 
stars (in the case of RHHW02 or LC03H). 

At this point, it should be noted that the aforemen- 
tioned yields are not the most appropriate for a discus- 
sion of the galactic ^°Fe/^^Al ratio. Indeed, the metallicity 
gradient observed in the Milky Way disk (-0.07 dex/kpc 
for oxygen and several other metals, see Hou et al. 2001 
and references therein) implies an average metallicity of 
around 2 Zq in the present-day disk. As already noted 
in several studies (e.g. Prantzos and Casse 1986), it is 
the yields of stars with such a metallicity that contribute 
mostly to the metal enrichment of the Milky Way today. 

Unfortunately, the works of RIIIIW02 and LC03 cover 
only solar metallicity stars, while the WW95 study con- 
siders a range of stellar metallicities below solar. One 
may, however, extrapolate from the trends obtained in 
the WW95 study at 2 Zq and scale accordingly the re- 
cent yields of RHHW02 and LC03. The WW95 yields for 
stars of initial metallicities Zq, 0.1 Zq and 0.01 Zq are 



displayed in Fig. 3, for ^^Al (upper panel) and for ^"Fe 
(lower panel), respectively. It can be easily see that the 
^"Fe yields are systematically proportional to the initial 
stellar metallicity for most stellar masses; the reason is 
that ^"Fe is mostly produced by neutron captures in the 
carbon shell and its yield is proportional to the initial ^^Fe 
amount. On the other hand, the yields of ^^Al are slightly 
higher at Zq than at 0.1 Zq. Part of ^^Al is produced in 
the H-shell by proton captures on initial ^^Mg and this 
does depend on initial metallicity; however, the bulk is 
produced in the C-shell (more than 80% in the 25 Mq 
star; see, e.g. Fig. 1 in Timmes et al. 1995) , where ^^Mg 
is produced by ^^C, itself resulting from the initial H and 
He of the star, and thus it is independent of the initial 
metallicity. 

One concludes then that at 2 Z© the ^°Fe yields of SNII 
(i.e. stars in the 12-25 Mq range) should be on average 
twice the corresponding ones at Zq, while the ^^Al yields 
should be only slightly higher than their counterparts at 
Zq. This implies in turn that the curves of ^°Fe/^^Al dis- 
played in Fig. 2 (corresponding to Zq stars) are in fact 
lower limits to the values expected from the galactic pop- 
ulation of massive stars This only exacerbates the dis- 
crepancy between the RHESSI result and the theoretical 
expectations from SNII, and makes the ^^Al contribution 
of WR stars even more important. Since the ^^Al yields 
of WR stars increase with metallicity approximately as 
Z^-^ or Z^(see below), they can easily match the increased 
^"Fe yields of SNII at 2 Zq and bring the average galactic 
^"Fe/^^Al ratio close to the RHESSI value. These qualita- 
tive considerations should be substantiated, of course, by 
self-consistent calculations of rotating stars at metallicities 
higher than Zq , extended as to cover all the advanced evo- 
lutionary phases, as well as the final explosion (see Heger 
et al. 2000 and Hirschi et al 2003 for preliminary results 
of such calculations). 

There is another way to understand the implications 
of the revised yields for the ^^Al sources, which does not 
involve ^°Fe. Indeed, observations of the Galactic 1.8 MeV 
line by different instruments converge to a value of 4 10^^ 
photons/cm^ /s , which corresponds to a steady state value 
of 2 Mq of ^^Al (produced per Myr) in the interstel- 
lar medium (e.g. Diehl et al. 1995, Prantzos and Diehl 
1996, Diehl and Timmes 1998). The average ^^Al yield in 
the recent calculations of SNII is 2.5 10~^ Mq (compared 
to 10~* Mq in WW95). Taking into account the average 
SNII frequency observed in Sbc galaxies like our own (~1- 
2 per century, Cappellaro et al. 2003) one sees that SNII 
may produce about 0.3-0.6 Mq of Al per Myr, i.e. about 
15-30 % of the total amount inferred from observations. 
The necessity for another source of ^^Al becomes obvious 
again. On the basis of non- rotating stellar models, Meynet 



In Timmes et al. (1995) a detailed model is presented for 
the Galactic disk, but the WW95 yields are used and no explicit 
mention is made of the galactic metallicity gradient. Thus, the 
ratio of 0.16 is obviously an underestimate, based on yields of 
stars with metallicity not exceeding Zq. 
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Fig. 4. Radial distributions of A126, star formation rate 
(SFR) and metallicity (Z) in the Milky Way disk. Upper 
panel: Data points with vertical error bars correspond to 
various tracers of the SFR, while the galactic metallic- 
ity profile of oxygen (with a gradient of dlog(O/H)=:-0.07 
dex/kpc) is shown by a solid line; the A126 profile, after 
an analysis of COMPTEL data by Knodlseder (1997), is 
shown (in relative units) by data points with vertical and 
horizontal error bars (the horizontal ones correspond to 
the adopted radial binning). Lower panel: If galactic A126 
originates mostly from WR stars, its radial distribution 
should scale with SFR * 7}-^ (points with vertical error 
bars, scaled from the upper panel), since the A126 yields of 
WR stars scale with Z^ '"^ at least (Vuissoz et al. 2003, cal- 
culations for rotating stars); however, the observed A126 
distribution (same points as in upper panel) is flatter than 
the expected one in that case. 

et al. (1997) show quantitatively that WR stars can indeed 
provide the bulk of galactic ^^Al. This is also supported 
by a different argument (Knodlseder 1999) concerning the 
similarity of the Galactic maps of Al and of ionizing pho- 
ton flux (provided only by the most massive stars, those 
that eventually become WR). Moreover, Knodlseder et al. 
(2001) point out that one of the prominent "hot-spots" of 
the COMPTEL 1.8 MeV map, the Cygnus region, is an 
association of very young massive stars, with no sign of 
recent supernova activity. 

Those arguments point towards WR stars as major 
sources of ^^Al in the Milky Way. However, the situation 
is far from being clear yet, because the WR stellar yields 
of ^^Al depend strongly on metallicity. In the case of non 
rotating stellar models that dependence is oc Z^, according 
to Meynet et al. (1997). The rotating models of WR stars, 
currently calculated by the Geneva group (Meynet and 
Maeder 2003) show that rotation considerably alleviates 
the need for high mass loss rates, while at the same time 
leading to the production of even larger ^^Al yields than 
the non-rotating models (Vuissoz et al. 2003); in that case, 
it is found that the Al yields of WR have a milder depen- 
dence on metallicity (oc Z^'^) than the non rotating ones. 



In both cases, that metallicity dependence of the ^^Al 
yields of WR stars, combined with the radial profiles of 
star formation rate (SFR) and of metallicity in the Milky 
Way (see Fig. 4, upper panel) suggest that the resulting 
radial profile of ^^Al should be much steeper than the one 
actually observed. The latter, derived from COMPTEL 
observations (Knodlseder 1997) appears in Fig. 4 (lower 
panel) and is clearly flatter than the product SFR*Z^-^ 
(as already noticed in Prantzos 2002). Similar conclusions 
are reached if the longitude, rather than radial, profiles of 
^^Al, metallicity and SFR are considered. 



3. Conclusion 

Contrary to a rather widely spread opinion, the recent 
RHESSI detection of radioactive ^^¥e in the Milky Way 
does not imply that ^^Al is mostly produced by supernova 
explosions. Recent theoretical results suggest that the ^"Fe 
line flux would then be close to the one of ^^Al (within 
a factor of two). Assuming that both the RHESSI results 
and the recent stellar nucleosynthesis results hold, another 
source of ^^Al should be found. 

Wolf-Rayet stars appear as natural candidates, in view 
of their absolute ^^Al yields (at least in the framework of 
the Geneva models: either with high mass loss rates and 
no rotation - Meynet et al. 1997 - or with mild mass loss 
rates and rotation - Vuissoz et al. 2003) and presumably 
low ^°Fe/^^Al ratios. However, the strong dependence of 
the ^^Al yields on metallicity suggests that the ^^Al emis- 
sivity should be steeply increasing in the inner Galaxy, 
while the COMPTEL observations clearly display a milder 
enhancement at small Galactic longitudes. 

Thus, almost twenty years after its discovery 
(Mahoney et al. 1982), the "^^PA emission of the Milky Way 
has not yet found a completely satisfactory explanation. 
Indeed, the recent observational (COMPTEL, RHESSI) 
and theoretical (RHHW02, LC03, Vuissoz et al. 2003) re- 
sults have made the puzzle even more complex than be- 
fore. The solution will obviously require progress in both 
directions. From the theory point of view, detailed nucle- 
osynthesis calculations of mass losing and rotating stars 
up to the final explosion in the mass range 12-100 M© 
and for metallicities up to 3 Zq will be required ; further- 
more, the uncertainties still affecting the reaction rates of 
^^Ne(a,n) (major neutron producer during He burning in 
massive stars) and ■^^Fe(n,7) will have to be substantially 
reduced. From the observational point of view, the radial 
distributions of both ^^Al and ^°Fe will be needed; such 
distributions will probably be available if the operation of 
ESA's INTEGRAL satellite is prolonged for a few years 
beyond its nominal 2-year operation. 
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